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Abstract 

The scattering function and radius of gyration of an ideal polymer network 
are calculated depending on the strength of the bonds that form the crosslinks. 
Our calculations are based on an exact theorem for the characteristic function 
of a polydisperse phantom network that allows for treating the crosslinks 
between pairs of randomly selected monomers as quenched variables without 
resorting to replica methods. From this new approach it is found that the 
scattering function of an ideal network obeys a master curve which depends 
on one single parameter x = (ak) 2 N/M, where ak is the product of the 
persistence length times the scattering wavevector, N the total number of 
monomers and M the crosslinks in the system. By varying the crosslinking 
potential from infinity (hard 5-constraints) to zero (free chain), we have also 
studied the crossover of the radius of gyration from the collapsed regime where 
R g ~ 0(1) to the extended regime R g ~ 0{yN). In the crossover regime the 
network size R g is found to be proportional to (N/M) 1 ' 4 . The latter result 
can be understood in terms of a simple Flory argument. 
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I. INTRODUCTION 



Sufficiently crosslinked macromolecules form solid-like rubber networks whose spectac- 
ular elastic properties are commonly believed to be of entropic origin. As a consequence 
polymeric networks are often successfully modeled as a set of independent random walks 
with the only restriction that the end points deform affinely under external stress. In a real 
network, however, permanent junctions between the macromolecules are randomly formed 
upon fabrication leading to a high degree of polydispersity. To meet with this complicated 
physical situation more sophisticated theoretical models are required. 

Recently various analytical studies have aimed at a statistical description of polymer 
networks taking randomness of the crosslinking positions into account. The mathematical 
challenge that arises in any of these approaches is that the permanent junction-points be- 
tween macromolecules are frozen and cannot be treated within the framework of Gibbsian 
statistical mechanics. This has first been realized by Edwards |1] and his replica formalism 
has by now become the standard approach in the field of polymer networks. Unfortunately 
the replica method forces strong approximations to remain analytically tractable. A recent 
work by Panyukov and Rabin || seems to overcome several difficulties and derives promising 
results for the elasticity. 

In this study we use a different route of thought. It has recently been shown |S| that 
for randomly crosslinked Gaussian structures substantial progress can be made by invoking 
quite different mathematical tools than replica field theory. The purpose of the paper is to 
report on further progress in this direction. Our working model is an ideal polymer chain 
of N monomers with M randomly selected pairs of monomers constrained to be in close 
neighborhood. The distance constraints are modelled by harmonic potentials. By varying 
the strength of the crosslinking potential we can continuously switch from a network situation 
with hard ^-constraints to the case of a free chain. 

Although the above model is highly idealized, the system is most interesting in its own 
right, since it is considered a first step towards a systematic theory of polymer networks 
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and gels. Moreover, the statistics of a huge macromolecule crosslinked to itself (randomly or 
not) has attracted a lot of attention recently because of its possible implications to protein 
structure reconstruction from NMR data P~pi|. Despite of the recent interest to the best of 
our knowledge these are the first exact results for ideal networks that have been reported. 



II. THE IDEAL NETWORK 

We adopt the minimal model of a huge Gaussian chain that is M times crosslinked to 
itself. In the Hamiltonian only terms that model chain connectedness and contributions due 
to crosslinking are retained. Complicating factors such as entanglements, excluded volume 
are deliberately neglected from the start. An appropriate Hamiltonian to begin with is 

= ± E(R, - r,-i) 2 + ± E(R e - %) 2 • (i) 

Za i=l Z£ e=l 

We have assumed N + 1 monomers whose locations in space are given by <i-dimensional 
vectors Rj (i = 0, 1, ...,N). Distance constraints exist between pairs of monomers labeled 
by i e and j e . For further use we introduce the inverse strength of the crosslinking potential 



e^ 2 



(2) 



as the mean squared distance between monomers that form the crosslinks measured in units 
of the persistence length a of the chain (figure 1). Limiting cases are given by z = (hard 
^-constraints) and z — > oo (free chain). The whole crosslinking topology is specified by a set 
of 2M integers C—{i e , j e }^_ v It has been shown || that the model in (0) is equivalent to the 
Deam-Edwards model without excluded-volume interaction if averages are understood in 
the following sense 

[utodRie-^.... 
(....)= lim J — r . (3) 

It is interesting to note that the entire range > z < oo can be treated with the same 
formalism. In || we have explicitly shown that all limits are mathematically well defined. 



Thus the here presented method is indeed able to cover the range of all crosslink strengths, 
i.e., from very soft to the quenched case. 

To model M uncorrelated crosslinks the distribution of frozen variables C is assumed to 
be uniform 



Other distributions are in principle possible but not considered in this investigation. As 
usual for systems with permanent constraints care must be taken in evaluating averages 
of physical quantities. The strategy here is not to start with the quenched average over 
the frozen variables by employing for instance the replica trick, but to keep explicitly all 
crosslink coordinates C during the calculation. Only at the very end the physical observable 
of interest is evaluated for a particular realization of C which is generated by the distribution 
in (f|). Clearly both approaches will give the same results if only self-averaging quantities 
are considered. 

The Hamiltonian in Eq. (1) together with the uniform distribution of crosslinks (4) 
defines our working model for the ideal network. 



In this section a brief review of the central mathematical theorem for the characteristic 
function of a Gaussian structure with internal (^-constraints is given, together with its ex- 
tension to arbitrary crosslinking potential z. The characteristic function for the problem is 
introduced as 



from which all expectation values can be obtained via differentiation. Simplifying notation 
has been adopted, where r,- = Rj — Rj_i (j = 1,...,N) denote bond vectors along the 
backbone of the chain, and E = (Ei, Ejv), r = (ri, ...,rjv) are ^-dimensional super- 
vectors with (i-dimensional vector components. Thus Zq(E; C) is also the partition function 




(4) 



III. CHARACTERISTIC FUNCTION 




(5) 
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of an ideal network in the presence of external fields Ej. Note that Zq(E; C) depends 
explicitly on all external fields contained in the vector E, as well as on all crosslink positions 
C. 

Without going into mathematical details, it is now possible to proof the following ana- 
lytically exact projection theorem || 

Z (E;C) = exp(-^Ei) , (6) 

where Ej_ is the length of the external field vector E projected perpendicular to the vector 
space spanned by "crosslink vectors" 

p 6 = (0,...,0, l,l,...,l,l ,0,...,0) . (7) 

i e + 1 to j e 

The above statement can be pictured in the following intuitive manner (figure 2). For each 
crosslink specified by the pair of integers 1 < i e < j e < N, form the corresponding N- 
dimensional vector p e , Eq. (^), where the l's are assumed to run from the (i e + l)th to 
the j e th position. The rest of the N components are filled with O's. The whole set of M 
vectors Pi, ...,Pm defines a characteristic vector space for the problem, say U. There is an 
unique decomposition of any field vector E parallel and perpendicular to U (figure 2), viz., 
E = E|| + Ej_. The operator that projects E on U can be constructed from p e for any 
realization of crosslinks C. It is given by VV + , where V is the N x M rectangular matrix 
associated with the crosslink vectors p e (e = 1, M), 

V = ( Pi ,...,pm) , (8) 



while V + is a generalized inverse of V |12[ . 



The projection theorem (^) is only valid for hard crosslinking constraints z = 0. A 
generalization to arbitrary crosslinking potential z is, however, possible using the methods 
of Ref. M. Only the final result for the characteristic function is quoted which reads 



2 (z, E; C) = exp 



>V E 2 ^ fx.E) 2 



2d V" § 1 + (z/w*) 



(9) 
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Here x e denotes any orthonormal basis associated with p e (e = 1, M), and w e are corre- 



sponding singular values fl2" |. For tetrafunctional crosslinks it was shown that w e is always 
positive. ^From (|9]) it is straightforward to obtain expressions for the radius of gyration and 
structure factor. 

By the scattering function (form factor) we mean density fluctuations [I3| normalized to 



one 



S (k, z; C) = (|p k 



(10) 



1 N 



N 2 A 

iV i,j=0 



where k is the scattering wavevector, and the average is over the measure in (j^) without 
having taken the z — > limit. In the following no distinction between N and N + 1 will 
be made since N is assumed to be large. From Eq. (|9D an exact expression for So can be 
obtained by introducing the external field vector E = fcCy, where 



(0. 



,0,1,1,. ..,1,1,0, 



0) 



(11) 



and k — Ikl. From 



i + 1 to j 
and (|T0D it is found 



So(Kz; C) 



N 



+ 



N 2 



^2 exp 



1 

N 

a 2 k 2 



(12) 



J ~ * 



\l + (z/w 2 e ), 



Similarly for the radius of gyration [13] 



'R g (z;Cy 



N 



N M 
i<j e =l 



[X. e Cij / 



(13) 



6 N 2 ^^l + (z/w 2 )- 

It is worthwhile to note that the applicability of these results is not limited to random 
networks since all crosslinking coordinates are still implicit in the formulas through x e and 
w e . By selecting different ensembles for C={i e ,j e }g£ 1 (random or not) any generalized 
Gaussian structure with internal crosslinking constrained can be treated by the same method. 
Further generalizations of the working model (^ to structures built from more than one chain 
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are feasible as long as the objects under investigation are simply connected. Otherwise the 
phantom character of the chains and the neglect of entanglements leads to derealization of 
those clusters of the network which are not connected via crosslinks. A similar percolation 



problem arises at the vulcanization transition [1^] which, however, is not an objective of this 
study. 

IV. DISCUSSION OF RESULTS 



In applying Eqs. (12) and (|13D to random networks, one has to deal with a sufficient 
number of monomers N and crosslinks M. Moreover, the positions of crosslinks C—{i e , je}e=i 
have to be chosen at random. The simplest scenario is to pick 2M integers i e ,j e (e = 
1,...,M) from the uniform distribution, Eq. (|]), defined on the interval [0, N]. For a 
given realization C we compute an orthonormal basis x e and singular values w e of p e (e = 



1, M). Any standard technique like singular value decomposition [K[ will suffice, for the 
orthonormalization process presents only a minor numerical task. We find that for number of 
monomers N > 10000 and crosslinks M > 200 fluctuations between different realizations of 
C differ by less than 1 percent. This also presents an estimate for the numerical uncertainty 
of the calculation. 

A. Scaling behavior of scattering function 

Within the framework of the Deam-Edwards model (z = 0) it was demonstrated || 
that the scattering function Sq of an ideal network is an universal function of wavevector k 
and mean crosslink density M/N as long as N and M are sufficiently large to ensure self- 
averaging. A scaling form for So is motivated by the following argument. For a linear polymer 
without crosslinks the scattering intensity S (x) depends only on the product x = k 2 R 2 , 
where R 2 = a 2 N/6 is the radius of gyration of the chain and Sq(x) = 2(e~ x — 1 + x)/x 2 
the Debye function [|T3 1 . In close analogy it was shown that for hard constraints (z = 0) 



the radius of gyration of an ideal network is given by R ~ 0.26 a N/M which suggests a 



scaling behavior, similar to that of linear chains without crosslinks. This scaling hypothesis 
for ideal networks was confirmed by our calculation based on the expressions in Eqs. ( |P2|) 
and (0). 

The numerically exact result for the scattering function is presented in figure 3. In the 
Kratky plot of figure 3, xS (x) has been evaluated as a function of y/x = kR s for z — 0. 
Independent of details of crosslinking topology C, all networks investigated fall on the same 
master curve (solid line). Statistical fluctuations between different networks were too small 
in the self-averaging regime to be seen on the scales used in figure 3. No indication of a 
power-law decay for intermediate wavevectors k or other simplifying feature was detected, 
besides the pronounced maximum in the Kratky plot at yfx ~ 2 which reflects the strong 
correlations of the monomers due to crosslinking. For comparison, the case of a linear 
polymer with z — > oo (Debye function, dashed line) was also computed from ([12]). 

The results of So for different values of crosslinking potential z are illustrated in figure 4 
for a network with N = 10000 and M = 200. By increasing the strength of the constraint 
from left to right, large deviations of So from ideal chain behavior (Debye function, left curve) 
arise on smaller and smaller length scales. For z = the network character persists down 
to even the shortest length scale k « 1 as a consequence of the high degree of crosslinking 
in the system (M/N = 0.02). For sufficiently large wavevectors all curves decay as k~ 2 as 
expected when the scanning wavelength becomes small compared to the mesh size of the 
network. Again, the master curve for So can be obtained by plotting the x-axis in units of 
R s . 

B. The collapse transition 

In figure 5 we have calculated the radius of gyration of a network of N monomers and 
M crosslinks as a function of y/z = e/a by use of Eq. (|I~3"|). From this investigation we can 
clearly distinguish three different scaling regimes 



S 



0.26 N/M, ife«ei 
R \ 2 

) - { 0.34 {e/ a) (N/M) 1/2 , if e 1 « s « e 2 (14) 



N/6 , if e » e 2 



with crossovers at e\ ~ ay N/M and £2 — ayMN . The plateau values in figure 5 correspond 
to the two extremes R 2 Ja 2 = 0.26 N/M (z -> 0) to the left and i^/a 2 = iV/6 (2 -> 00) to 
the right. 

In particular our investigation showed that the cases z = (hard constraints) and z — 1 
(constraints of the order of the persistence length a) only differ by a numerical prefactor 
which varies from 0.26 for z = to about 0.27 for z = 1. From this we conclude that an 
ideal network subject to uncorrelated crosslinking constraints is collapsed in a sense that 
its size is proportional to the square root of N/M. Therefore Rja — 0(1) since M/N is 
the mean crosslink density in the network and of order unity in the thermodynamic limit 
N,M — > 00. This finding seems to be at variance with current speculations regarding the 
collapse transition of macromolecules []10[, where it was argued that a critical number of 



crosslinks M > M c ~ N/ log iV will force the system to collapse. Our result for z = is in 
agreement with recent Monte Carlo simulations by Kantor and Kardar [ I I ] who found for 
the mean squared end-to-end distance R 2 /a 2 ~ 1.5 N/M. This suggests the same one to six 
ratio for (R g /R) 2 in ideal networks as for linear polymers without crosslinking constraints 
and excluded- volume interaction [13]. We believe that the discrepancy between our exact 
results and these suggested in [HJ are due to the approximations used there. 



Conversely, a free chain (z — > 00) is an extended object with R g /a ~ 0{yN). Between 
the collapsed and the extended regime we find a smooth crossover with (R g /a) 2 being pro- 



portional to {e/a)<\jN/M. Remarkably this is the same scaling as for randomly branched 
polymers without excluded volume interaction. In the following we discuss the different 
conformational states of the system in terms of simple scaling arguments. 
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C. Flory estimates 



For completeness we first note that the free chain regime is trivial since R 2 /a 2 = N/6 is 
an exact solution of (13) for z — > oo. To understand the scaling behavior of R s in the other 
two regimes the free energy of the Hamiltonian in Eq. (1) is discussed within Flory theory. 
The connectivity term in (1) models the standard entropic elasticity of a Gaussian chain, 
i.e., R 2 /(Na 2 ) + Na 2 /R 2 , where R is a measure of the size of the system. The first term 
accounts for stretching, whereas the second term describes the response due to compression 

An estimate of the crosslink term in (1) requires more attention. First we consider soft 
crosslinks when e » a. In this regime the second term of the Hamiltonian is estimated by 
M(R/e) 2 , because the mean squared distance between a pair of constrained monomers is of 
order e 2 . The relevant part of the total Flory free energy then is 

No 2 MR 2 

+ — ■ (15) 

Minimization of the free energy with respect to R yields the scaling relation R/a ~ 
(e/a) 1 ' 2 (N/M) 1 / 4 in agreement with (14). The appearence of the branched polymer ex- 



ponent 1/4 can be assigned to the change of connectivity of the chain when e ~ 0(\jNM). 

The case of hard crosslinks e ~ 0(a) is more difficult to obtain. We picture the system as 
a coarse-grained random walk over the M crosslinks with an effective step length proportional 
to N/M, i.e., the mean number of monomers between crosslinks. From this mean-field 
argument the crosslink term is estimated to be of the order M[R 2 /(a 2 N/M)}. The latter 
expression has the effect that it tries to shrink the chain upon cost of confinement entropy. 
A suitable Flory free energy is given by 

Na 2 M 2 R 2 

^°~l*r + W (16) 

^From there (R/a) ~ (N/M) 1 / 2 as was shown exactly for the network regime. 
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FIGURES 

FIG. 1. Modelling of a soft crosslink between two monomers i e and j e . A spring constant e — > 
leads to hard (5-constraints, e — > oo to free chain behavior. 

FIG. 2. External field vector E projected parallel and perpendicular to U. The vector space U 
is spanned by the vectors p e defined in Eq. (7). The sketch is for two crosslinks e = 1,2. 

FIG. 3. Kratky plot of xSq(x) for the ideal network (z = 0, solid line) and linear chain [z — > oo, 
Debye function, dashed line) as a function of y/x = fcl? g in d = 3. The open circles show a second 
realization of C to demonstrate self-averaging. Note that we measure length scales in terms of R g 
and use a normalization So(0) = 1 different to that in the experimental literature. 

FIG. 4. Structure function Sq for different crosslinking potentials z. From right to left 
z = 0, 10, 10 2 , 10 3 , 10 4 , 10 5 , oo for a network size of N = 10000 and M = 200 plotted in dimension- 
less units q = ka/y/2d. The left curve (z — ► oo) is the Debye function of a linear chain. 

FIG. 5. Radius of gyration of a polymer chain as a function of \fz = e/a. For the three solid 
curves the number of monomers was varied from top to bottom iV = 5000, 10000, 20000; M was 
kept constant at 200. The short dashed line shows a network with iV = 20000 and M = 400. 
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